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Abstract. We study the fate of static fundamental charges in the thermodynamic limit 
from Monte-Carlo simulations of SU{2) with suitable boundary conditions. 



1. Introduction 

In QED, the charge of a particle is of long-range nature. It can exist 
because the photon is massless. Localized objects are neutral like atoms. 
Within the language of local field-systems one derives more generally that 
every gauge-invariant localized state is singlet under the unbroken charges 
of global gauge invariance. Thus, without (electric) Higgs mechanism, QED 
and QCD have in common that any localized physical state must be charge- 
less/colorless. 

The extension to all physical states is possible only with a mass gap. 
Without that, in QED, non-local charged states which are gauge-invariant 
can arise as limits of local ones which are not. The Hilbert space decomposes 
into the so-called superselection sectors of the physical states with different 
charges. With a mass gap in QCD, on the other hand, color-electric charge 
superselection sectors cannot arise: every gauge-invariant state can be ap- 
proximated by gauge-invariant localized ones (which are colorless). One 
concludes that every gauge- invariant state must also be a color singlet. 

On the other hand, charged states are always possible with suitable 
boundary conditions in a finite volume. This allows to study their fate in 
the thermodynamic limit from Monte-Carlo simulations on finite lattices. In 
an Abelian theory for example, anti-periodic (spatial) boundary conditions 
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can be used to force the system into a charged sector in the infinite volume 
Umit Q. The (Higgs vs. Coulomb) phases of the non-compact Abelian Higgs 
model can be distinguished in this way. And by duality, via the Z gauge 
theory, the magnetic sectors of compact U{1) follow an analogous pattern. 
The difference in free energy of the anti-periodic vs. the periodic ensemble 
thereby tends to zero or a finite value for the (magnetic) Higgs or Coulomb 
phases, respectively. 

In pure SU{N) gauge theory, one expects the free energy Fq{T,L) of 
a static fundamental charge in a 1/T x box, for L— >oo, to jump from 
oo to a finite value at T = Tc reflecting the deconfinement transition. The 
Polyakov loop P is commonly used to demonstrate this in lattice studies. 
If (P) = e"^'/"^, the center symmetric (broken) phase gives for Fq an infi- 
nite (finite) value. However, the periodic boundary conditions (b.c.) within 
which (P) is measured are incompatible with the presence of a single charge 
also in this case. And, like any Wilson loop, (P) is subject to UV-divergent 
perimeter terms, such that (P) = at all T as the lattice spacing a— >0. 

Following it is possible, however, to measure the gauge- invariant, 
UV-regular free energy of a static fundamental charge ^ , and show that 
it has the expected behaviour, dual to that of temporal center flux ||5|. The 
preparation of suitable b.c.'s to achieve this is a little indirect. 

2. Twist vs. Electric Flux Sectors in SU{2) 

For the different sectors relevant to the confinement transition in pure 
SU{N) gauge theory, one needs to distinguish between the finite volume 
partition functions of two types. 

First, 't Hooft's twisted boundary conditions fix the total number of 
Zjv-vortices modulo N that pierce planes of a given orientation. On the 
4-dimensional torus there are A^^ different such sectors correspond- 
ing to the 6 possible orientations for the planes of the twists. Without 
fields that faithfully represent the center Ztv of SU{N), the structure is 
G = SU{N)/Zn with first homotopy 7ri(G) = Zat. The A^^ inequiv- 
alent choices for imposing (twisted) boundary conditions on the gauge 
potentials A therefore correspond to the classification of the bundles, by 
their Zjv-vortex numbers, according to the harmonic 2- forms over with 
TTi{G) = Zat coefficients, the 2nd de Rahm cohomology group H'^(T'^, I-n)- 

At finite temperature T > the A^^ possible twists come in two classes: 
3 temporal ones classified by a vector k G Z^r, and 3 magnetic ones by 
fn e Z^r, see Fig. 0. Magnetic twist is defined in purely spatial planes and 
fixes the conserved, Zjv-valued and gauge-invariant magnetic flux rh in the 
perpendicular directions. 
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3d-line defect: Z2-vortox, maps circle to a non-contractible 
loop in IRP(2), same happens in 50(3) = 5(7(2)/Z2. 



Figure 1. 3d Z2-vortex as in, e.g., nematic liquid crystals (left). Assume the 4d-vortices 
can lower their free energy by spreading (right): they can do that in the 3 spatial m-planes 
at all T (top), while they are squeezed in the 3 temporal fc-planes at large T (bottom). 

The different clioices of twisted b.c.'s lead to sectors of fractional Chern- 
Simons number [v + k ■ fh/N) Q with states labelled by m, i^), where 
G Z is the usual instanton winding number. These sectors are connected 
by homotopically non-trivial gauge transformations r2[A;,i/], 

9.[k' ,v']\k,m,u) = \k+k' ,m,v+v') . (1) 

A Fourier transform of the twist sectors Zk{k,m,v)., which generalizes the 
construction of 0-vacua as Bloch waves from z^-vacua in two ways, by re- 
placing V ^ [v + k ■ fh/N) for fractional winding numbers and with an 
additional Z^-Fourier transform w.r.t. the temporal twist k, yields, 

Ze{e,m,e) = e-^^=(^^'-'^) = ^ E e"^"^''") m, z.) . (2) 

k, V 

Up to a geometric phase Lo{k, u) = lire - k/N + 9{v+k ■ fh/N), the states in 
the new sectors are then invariant under the non-trivial Q.[k, u] also, 

VL[k,u\\e,fh,9) = eyip{iuj{k,v)}\e,fh,9) . (3) 

Their partition functions Zf, are classified, in addition to their magnetic 
flux fh and vacuum angle 9, by their Z^v-valued gauge-invariant electric 
flux in the e-direction |Q. Here, we do not consider finite 9 which we omit 
henceforth. Recall the following points for details of which we refer to Q: 





Figure 2. Cubes with one, two, and three LxL planes dual to the stacks of plaquettes 
flipped for temporal twist k = (0, 0, 1), k = (0, 1, 1), and k = (1, 1, 1), from left to right. 



(i) The twisted partition functions, relative to the periodic ensemble 
Zfc(0, 0), are expectation values of combinations of 't Hooft loops W of 
maximal size in (/i, z^)-planes dual to the planes of the twists, 

Zk{k,rh)/Zk{m = {Wj^%) ■ (4) 

In particular, the temporal A;-twists correspond to expectation values of 
spatial 't Hooft loops. The Z^-Fourier transform of Eq. (^) exhibits their 
Kramers- Wannier duality with the electric flux sectors which are expecta- 
tion values of Polyakov loops in the no- flux ensemble Ze(0, 0) (see below). 

(ii) Note also that the no-flux ensemble in a finite volume is manifestly 
different from the periodic ensemble, e.g., for SU{2) one has, 

ZeiO, 0) = (l + 3(Tyi--t ) + m^t^.t ) + (W^sTat )) (5) 

The combinations of spatial 't Hooft loops needed to compute this, or any 
of the electric flux partition functions Ze{e,0), are sketched in Fig. ^. 

(iii) From the gauge-invariant definition of the Polyakov loop P{x) in 
presence of temporal twist ||7|, it is relatively simple but important to verify 
that the electric-flux partition functions are indeed expectation values of 
P's in the no-flux ensemble which follow the general pattern, 

Ze(o,o) » "i.r |i, (H';»»;,)-o, r>r, ^' 

with e G again. For ej^O this therefore yields the free energy Fq{T, L) = 
Ff>{e, 0) of one static fundamental charge in the volume with b.c.'s such 
that its electric flux is directed towards its 'mirror' (anti)charge in the 
adjacent volume along the direction of e. Also note that the operator in the 
expectation value of has no perimeter, is UV-regular, and one can see 
in Fig. ^ that there is no Coulomb term for small volumes either. 

Of course, Eq.(|6|) reflects the different realizations of the electric Z^ 
center symmetry in the respective phases. As compared to spin correlations 
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Figure 3. Partition function and free energy of temporal twist (top) and electric flux 
(bottom), both over temperature (left) and finite size scaling variable x (right). 



of the form {sg s^j^Lg.) in the 3d-Ising model with interfaces, the Polyakov 
loops in ^ are the corresponding variables in SU{2), whose behavior as a 
function of temperature is reversed. 

The dual of the Ising model on the other hand is the 3d Z2 gauge theory. 
Interfaces in the first are Wilson loops in the latter. Through duality, the 
expectation value of a Z2-Wilson loop along C is expressed as ratio of 
Ising-model partition functions with and without antiferromagnetic bonds 
at links dual to a surface 0, the Z2-interface, bounded by C, e.g., see 

In SU (2), the objects dual to the Polyakov loop correlations, or the electric 
fluxes in (|6|), are spatial 't Hooft loops. Via universality, these are the Z2- 
Wilson loop analogues. And their expectation values are calculated on the 
lattice in much the same way, by flipping a coclosed set 0* of plaquettes 
dual to some surface subtended by the loop C, 

{WsumiC))= j,Ue.p-(3j:^TrU^ ' = j /5 , □ ' 

In both cases the surface is arbitrary except for its boundary C = dO,. 
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A spatial 't Hooft loop of maximal size Lx L, living in, say, the (x, y) plane 
of the dual lattice, is equivalent to an odd number of flipped plaquettes in 
every {z, t) plane of the original lattice. This enforces twisted b.c.'s in {z, t). 
Combining such loops yields the other fc-twist sectors, c./.. Fig. |2|. 

The temperature dependences of the partition functions for temporal 
twist k = (0, 0, 1) and electric flux e = (0, 0, 1) as calculated in |^] are com- 
pared in Fig. p. Their dual behavior is obvious: for L^cc, both approach 
step functions jumping from 1 to and to 1, respectively, as Tc is crossed 
(from below). Near the phase transition, this behavior is determined by 
critical exponents and likewise universal amplitude ratios of the 3d-Ising 
class. For the larger spatial lattice sizes, the fits in the left part of Fig. |3| 
might look rather daring at first. However, each of the two families of curves 
shown there really represent one of the unique functions in the right part 
which fit all the data. This is possible due to finite size scaling. 



3. Finite-Size Scaling 



Finite size scaling (FSS) laws are based on the observation that the length 
of the system L and the correlation lengths ^ that diverge in the thermo- 
dynamic limit are the only relevant length scales in the neighborhood of 
the transition. In particular, as the critical point is approached, the finite 
lattice spacing a becomes less and less important. For the continuous (2nd 
order) transition of SU{2) in the 3d-Ising class, with £,±{t) = ^±\t\'^ for 
reduced temperature i ^ 0, we use t = T/Tc — 1 and u = 0.63. And as for 
the ratios of 3d-Ising model partition functions Za/Zp with anti-periodic vs. 
periodic b.c.'s |^], we assume the L, T dependence of the various temporal 
twist sectors, denoting Zk{i) = (l^™pat ) for i = 1, ..3 orthogonal fc-twists, 
to be governed by simple FSS laws, 

Zk{i) = /W(x) , with X = iLTcltl^ cx L/C±{t) , for t ^ . (9) 

We then observe that our results over x for all different lattice sizes nicely 
collapse on a single curve, c./.. Fig. ^. In the high temperature phase above 
Tc, the large-x behavior, — ln/'''(x) = ctq^ x'^ + . . ., refiects the dual area 
law for (i = 1, ..3) large spatial 't Hooft loops. Their dual string tension, 

is the Ising analogue of the interface tension (yi{t) = R/^'^{t) ~ \t\'^'^ below 



Tc, where R ~ 0.104 is a universal ratio |1C , [Til . In addition, the universality 
hypothesis relates the ratio of the correlation lengths for the Polyakov loops 
in SU{2) to the correlation lengths of the spins in the 3d-Ising model, as 



measured in [11|, 



^SU{2) i^SU{2) _!_ Rising ^^ising r;^ ]^ gg (]^]_) 
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Figure A. Free energy -ffc(l) = — InZfc(l) of one temporal twist in 5*17(2), same data as 
in Fig. pi compared to that of an interface in the Ising model -Ffe(l) = — \n{Za/Zp) with 
a;ising ~ — 1.88(2)a;sj7(2) (left); 1-3 orthogonal fc-twists in SU{2) (T>Tc) vs. Ising with 
1-3 anti-periodic directions {T<Tc), with square- root ratios in the fits (right), from llSl. 



Together with the large-x asymptotics ofi = 1, ..3 electric fluxes, —\nZe{i) = 
(7q'|x| + . . ., below Tc, this relates the string tension amplitude fio = ctq^ 
below Tc to its dual counterpart do = (Tq^' above Tc, 

a{T) = a^TcW T , with ao = ^ao/R+ , i?+ = i^eVCl - 0.4 . (12) 

Within the accuracy of our results, the SU (2) data is fully consistent with 
these ratios. The quite impressive range of universal behavior can be appre- 
ciated in comparing the SU (2) temporal twist with the interface free energy 
— \\iZa/Zp in the 3d- Ising model [^], see Fig. ^. Once a non-universal 
constant of proportionality relating both FSS variables is fixed, the scaling 
functions appear to be identical for the whole range of the SU{2) data 
(with high and low temperature phases interchanged). The same agreement 
is observed for all f^^\x), i.e., also between 2 and 3 orthogonal twists in 
SU{2) and Ising model with anti-periodic b.c.'s in 2 and 3 directions. 

4. String Formation 

In the low temperature phase, the formation of electric flux strings is ex- 
pected. The signature for this are square-root ratios 1 : \/2 : \/3 of the string 
tension amplitudes (Tq' for i = 1,2 and 3 orthogonal fluxes. At T = the 
singificance for such a behavior, as compared to the ratios 1:2:3 expected 
for isotropic fluxes, was assessed in the pioneering study of Ref. [ |12| . 

Above Tc on the other hand, the same square-root ratios for the dual 
string tension amplitudes, 

: erf -.aj^^ ~ 1:V2:V3, (13) 
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signal the formation of interfaces with minimal area. These ratios are well 
confirmed for spatial 't Hooft loops in orthogonal planes within the accuracy 
of our SU{2) data Q, and with considerably higher accuracy also for the 
3d-Ising model (for T < Tc) with anti-periodic b.c.'s in 1,2 and 3 directions 
12], which are also compared in Fig. El 



5. Conclusions 



To summarize, we have studied the finite volume partition functions in the 
sectors of pure SU{2) of two types: Using 't Hooft's twisted boundary con- 
ditions we first measured the free energies of ensembles with odd numbers 
of Z2 center-vortices through temporal planes. From combinations of these 
we then obtained the sectors of gauge-invariant electric flux, and demon- 
strated explicitly that, below Tc, their free energy diverges linearly with 
the length L of the system. Because spatial twists share their qualitative 
low-T behavior with the temporal ones considered so far, the free energy 
of the magnetic fluxes must vanish just as that of temporal twist. This is 
the magnetic Higgs phase with electric confinement of pure SU{2). 

At criticality all free energies rapidly approach their finite L ^ 00 limits. 
Above Tc, the electric-fiux free energies vanish in the thermodynamic limit. 
The transition is well described by simple finite size scaling laws of the 3d- 
Ising class. The ratios of the (dual) string tension amplitudes for 1,2 and 3 
orthogonal (spatial 't Hooft loops) electric fiuxes (above) below Tc indicate 
the formation of diagonal (interfaces) flux strings. 

Both of us would like to express our warm thanks to S. Olejnik and 
J. Greensite for their great job organizing this stimulating workshop. 
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